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ABSTRACT 

We present new accurate near-infrared (NIR) spheroid (bulge) structural parameters obtained by two- 
dimensional image analysis for all galaxies with a direct black hole (BH) mass determination. As expected, 
NIR bulge luminosities Lbui and BH masses are tightly correlated, and if we consider only those galaxies with se- 
cure BH mass measurement and accurate Lbui (27 objects), the spread of MBH-ibui is similar to MBH-Ce, where ag 
is the effective stellar velocity dispersion. We find an intrinsic rms scatter of ~ 0.3 dex in logMeH- By combining 
the bulge effective radii R^, measured in our analysis with <Je, we find a tight linear correlation {rms ~ 0.25 dex) 
between Mbh and the virial bulge mass (oc Recrl), with (MBH/Mbui) ^ 0.002. A partial correlation analysis shows 
that Mbh depends on both and Re, and that both variables are necessary to drive the correlations between Mbh 
and other bulge properties. 

Subject headings: black hole physics - galaxies: bulges - galaxies: nuclei - galaxies: fundamental parameters 



1. INTRODUCTION 

Central massive black holes (BHs) are now thought to reside 
in virtually all galaxies with a hot spheroidal stellar compo- 
nent (hereafter bulge). Such BHs seem to be a relic of past 
quasar activity (e.g., Soltan 1982; Marconi & Salvati 2002; Yu 
& Tremaine 2002; Aller & Richstone 2002) and related to host 
galaxy properties, with the implication that BH and galaxy for- 
mation processes are closely linked. Previous work has shown 
that BH mass Mbh is correlated with both blue luminosity LB.bui 
and bulge mass Mbui, although with considerable intrinsic scat- 
ter {rms ^ 0.5 in logMBn; Kormendy & Richstone 1995). How- 
ever, Mbh and the bulge effective stellar velocity dispersion ag 
correlate more tightly {rms ^ 0.3) than Mbh -i-B.bui (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000). The smaller scatter of the 
Mbh-o^c correlation suggests that the bulge dynamics (or mass), 
rather than luminosity, is the agent of the correlation. But the 
smaller spread relative to MBH-^bui appears to be an artefact 
of the manipulations necessary to derive Lbui- Indeed, recent 
work has shown that when bulge parameters are measured with 
more accuracy [e.g. profile fitting rather than average correc- 
tion for disk light (Simien & de Vaucouleurs 1986)], the result- 
ing scatter is comparable to that of Mbh-Cp (McLure & Dun- 
lop 2002; Erwin et al. 2003). The correlation between Mbh 
and bulge light concentration also has a comparably low scatter 
(Graham et al. 2001). Nevertheless, there are strong indications 
that LB.bui of the brightest elliptical galaxies, for which decom- 
position issues are unimportant, deviate significantly from the 
A^BH-^bui relation (Ferrarese 2002). Hence, longer wavelengths 
may also be necessary to better define the intrinsic scatter in 
MBH-Lbui compared to that of Mbh-C(.. 

In this paper, we reexamine the MBH-Lbui correlation by ac- 
curately measuring the bulge luminosity in the near-infrared 
(NIR) for all galaxies with a well-determined Mbh- All previ- 
ous studies have used optical light (B or R) to test the MBH-Lbui 
relation, but NIR light provides a clear advantage over the op- 
tical: it is a better tracer of stellar mass and less subject to the 
effects of extinction. If the physical correlation is between the 
BH mass and bulge mass, the NIR correlations MBH-Lbui should 
be tighter than those in the optical, because of the smaller vari- 



ation of M/L ratio T with mass (e.g., Gavazzi 1993). More- 
over, we use a two-dimensional (2D) bulge/disk decomposition 
to determine bulge parameters, an improvement on earlier work 
which applied ID fits only. Here we construct the largest pos- 
sible sample, by considering all galaxies which have been used 
for the MBH-Cf and MBH-i-B.bui correlations. In §2 we present 
the sample of galaxies with direct dynamical BH mass mea- 
surements, and in §3 describe the images and the 2D bulge/disk 
decomposition applied to them. Finally, in §4 we discuss the 
results of the analysis. 

2. THE SAMPLE 

To date, there are 37 galaxies with direct gas kinematical or 
stellar dynamical determination of the central BH mass. These 
galaxies have been compiled and made into a uniform sam- 
ple (e.g., for distances) by a number of authors (e.g., Mer- 
ritt & Ferrarese 2002; Tremaine et al. 2002; hereafter MF02 
and T02, respectively). We adopt the data from the recent pa- 
per by T02 with some modifications and additions. The data 
in Columns 1-5 and 9 of Table 1 are from the compilation 
by T02 and the reader can refer to that paper for more de- 
tails. Differently from T02, when galaxy distances from sur- 
face brightness fluctuations (Tonry et al. 2000) are not avail- 
able, we use recession velocities corrected for Virgocentric in- 
fall from the LEDA database (http://leda.univ-lyonl.fr/) with 
Ho = 70 km S-' Mpc"' . In a few cases, we also consider BH 
mass estimates from different papers than those used by T02; 
thus, in Col. 6, we indicate the appropriate references. With 
respect to the 3 1 galaxies considered by T02 we add: Cygnus 
A (Tadhunter et al. 2003), M81 (Devereux et al. 2003), M84 
(Bower et al. 1998), NGC 4594 (Kormendy 1988), Centaurus 
A (Marconi et al. 2001) and NGC 5252 (Capetti et al. 2003). 

Following MF02, we divide the galaxies into two groups. In 
the first group, we place all the galaxies which have a secure 
BH mass measurement and an accurate determination of the 
bulge NIR luminosity. We consider 'secure' those BH masses 
for which the black hole sphere of influence, R^h = GMbh/ct^ 
(column 7 of Table 1), has been clearly resolved, i.e., Nres = 
IR^a/Rres > L where Rres is the spatial resolution of the obser- 
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Table 1 

Galaxy Sample. 



Galaxy Type D a^, Mbh(+,-) Ref «bh N,es Mb M, Mh Mk Re Mw 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 



Group 1 



NGC4258 


Sbc 


7.2 


130 


3.9(0.1,0.1) 


X 


10' 


m-1 


0.28 


71 


-17.2 


-20.9 


-22.0 


-22.4 


0.92 ±0.23 


1.1±0.3 


X 10"* 


M87 


EO 


16.1 


375 


3.4(1.0,1.0) 


X 


10" 


g-2 


1.33 


33 


-21.5 


-24.6 


-25.2 


-25.6 


6.4 ±1.6 


6.2 ± 1.7 


X 10" 


NGC3n5 


SO 


9.7 


230 


9.1(9.9,2.8) 


X 


10* 


s-3 


1.57 


15 


-20.2 


-23.5 


-24.2 


-24.4 


4.7 ±1.2 


1.7 ±0.5 


X 10" 


NGC4649 


El 


16.8 


385 


2.0(0.4,0.6) 


X 


10" 


s-4 


0.71 


14 


-21.3 


-24.9 


-25.5 


-25.8 


8.1±2.0 


8.4 ±2.2 


X 10" 


M81 


Sb 


3.9 


165 


7.6(2.2, 1.1) 


X 


10' 


g-5 


0.63 


13 


-18.2 


-23.1 


-23.9 


-24.1 


3.4 ±0.9 


6.4 ± 1.8 


X 10'° 


M84 


El 


18.4 


296 


1.0(2.0,0.6) 


X 


10" 


g-6 


0.55 


11 


-21.4 


-24.7 


-25.8 


-25.7 


8.2±2.1 


5.0 ± 1.4 


X 10" 


M32 


E2 


0.8 


75 


2.5(0.5,0.5) 


X 


10^ 


s-7 


0.49 


9.7 


-15.8 


-18.9 


-19.7 


-19.8 


0.24 ±0.06 


9.6 ± 2.6 X 10* 


CenA 


SO 


4.2 


150 


2.4(3.6,1.7) 


X 


10* 


g-8 


2.25 


9.0 


-20.8 


-23.8 


-24.3 


-24.5 


3.6 ±0.9 


5.6 ±1.5 X 10'° 


NGC4697 


E4 


11.7 


177 


1.7(0.2,0.1) 


X 


10* 


s-4 


0.41 


8.2 


-20.2 


-23.9 


-24.5 


-24.6 


9.1 ±2.3 


2.0 ±0.5 X 10" 


IC1459 


E3 


29.2 


340 


1.5(1.0,1.0) 


X 


10' 


s-9 


0.39 


7.8 


-21.4 


-24.8 


-25.3 


-25.9 


8.2 ±2.0 


6.6±1.8 


X 10" 


NGC5252 


SO 


96.8 


190 


1.0(0.2,0.4) 


X 


10' 


g-10 


0.25 


5.1 


-20.8 


-24.4 


-25.2 


-25.6 


9.7 ±2.4 


2.4 ±0.9 X 10" 


NGC2787 


SBO 


7.5 


140 


4.1(0.4,0.5) 


X 


10' 


g-11 


0.25 


5.0 


-17.3 


-20.4 


-21.1 


-21.3 


0.32 ±0.08 


4.4 ± 1.2 X 10' 


IN KJK^^■Jy^ 


oa 


O fi 


240 


1 .Ul,i.U, u. /, 


X 


in' 


s-12 


I.J / 


^ n 

D.U 


-Zi. J 


-24 2 




0^ A 


^ 1 _|_ 1 


2.0 ±0.5 X 10" 


NGC3608 


E2 


22.9 


182 


1.9(1.0,0.6) 


X 


10* 


s-4 


0.22 


4.4 


-19.9 


-23.4 


-24.0 


-24.1 


4.3±1.1 


9.9 ±2.7 X 10"* 


NGC3245 


SO 


20.9 


205 


2.1(0.5,0.5) 


X 


10* 


g-13 


0.21 


4.2 


-19.6 


-22.4 


-23.1 


-23.3 


1.3 ±0.3 


3.9 ± 1.0 X 10"* 


NGC4291 


E2 


26.2 


242 


3.1(0.8,2.3) 


X 


10* 


s-4 


0.18 


3.6 


-19.6 


-23.1 


-23.8 


-23.9 


2.3 ±0.6 


9.5 ±2.5 X 10"* 


NGC3377 


E5 


11.2 


145 


1.0(0.9,0.1) 


X 


10* 


s-4 


0.38 


3.6 


-19.0 


-22.7 


-23.5 


-23.6 


5.4 ±1.3 


7.8±2.1 


X 10«* 


NGC4473 


E5 


15.7 


190 


1.1 (0.4,0.8) 


X 


10* 


s-4 


0.17 


3.4 


-19.9 


-23.1 


-23.6 


-23.8 


2.8 ±0.7 


6.9 ± 1.9 


X 10'° 


CygnusA 


E 


240 


270 


2.9(0.7,0.7) 


X 


10" 


g-14 


0.15 


2.9 


-21.9 


-26.4 


-26.9 


-27.3 


31 ±8 


1.6± 1.1 


X 10'^ 


NGC4261 


E2 


31.6 


315 


5.2(1.0, 1.1) 


X 


10* 


g-15 


0.15 


2.9 


-21.1 


-24.6 


-25.4 


-25.6 


6.5 ± 1.6 


4.5 ± 1.2 


X 10" 


NGC4564 


E3 


15.0 


162 


5.6(0.3,0.8) 


X 


10' 


s-4 


0.13 


2.5 


-18.9 


-22.5 


-23.3 


-23.4 


3.0 ±0.7 


5.4± 1.5 


X 10"* 


NGC4742 


E4 


15.5 


90 


1.4(0.4,0.5) 


X 


10' 


s-16 


0.10 


2.0 


-18.9 


-22.1 


-22.8 


-23.0 


2.0 ±0.5 


1.1 ±0.3 


X 10"* 


NGC3379 


El 


10.6 


206 


1.0(0.6,0.5) 


X 


10* 


s-17 


0.20 


1.9 


-19.9 


-23.1 


-23.7 


-24.2 


2.9 ±0.7 


8.5 ±2.3 


X 10"* 


NGC1023 


SBO 


114 


205 


4.4(0.5,0.5) 


X 


10' 


s-18 


0.08 


1.6 


-18.4 


-22.6 


-23.3 


-23.5 


1.2 ±0.3 


3.4 ±0.9 


X 10'° 


NGC5845 


E3 


25.9 


234 


2.4(0.4,1.4) 


X 


10* 


s-4 


0.15 


1.4 


-18.7 


-22.0 


-22.7 


-23.0 


0.50 ±0.12 


1.9 ±0.5 


X 10'° 


NGC3384 


SO 


11.6 


143 


1.6(0.1,0.2) 


X 


10' 


s-4 


0.06 


1.2 


-19.0 


-21.7 


-22.3 


-22.6 


0.49 ±0.12 


7.0 ± 1.9 


X 10» 


NGC6251 


E2 


107.0 


290 


6.1(2.0,2.1) 


X 


10* 


g-19 


0.06 


1.2 


-21.5 


-25.4 


-26.4 


-26.6 


11±3 


6.7 ±1.8 


X 10" 


Group 2* 


MilkyWay 


SBbc 


0.008 


103 


4.1(0.6,0.6) 


X W' 


s-20 


42.9 


1714 


-17.6 


-22.0 


-22.2 


-22.3 


0.70 ±0.20 


5.2 ±2.5 


X 10» 


M31 


Sb 


0.8 


160 


4.5(4.0,2.5) 


X 10' 


s-21 


2.05 


41 


-19.0 


-21.8 


-22.5 


-22.8 


1.0±0.3 


1.9 ±0.5 


X 10"* 


NGC1068 


Sb 


15.0 


151 


8.3(0.3,0.3) 


X 10''' 


m-22 


0.02 


2.7 


-18.8 


-23.6 


-24.7 


-25.0 


3.1 ±0.8 


5.0 ± 1.4 


X 10"* 


NGC4459 


SO 


16.1 


186 


7.0(1.3, 1.3) 


X 


10' 


g-11 


0.11 


2.2 


-19.1 


-23.9 


-24.2 


-24.5 


15±4 


3.6 ± 1.0 


X 10" 


NGC4596 


SBO 


27.9 


152 


7.8(4.2,3.3) 


X 


10' 


g-11 


0.11 


2.1 


-20.6 


-23.0 


-23.7 


-23.8 


1.6 ±0.4 


2.6 ± 0.7 


X 10'° 


NGC7457 


SO 


13.2 


67 


3.5(1.1, 1.4) 


X 


10" 


s-4 


0.05 


1.0 


-17.7 


-21.3 


-22.0 


-21.8 


4.8 ±3.5 


1.5± 1.1 


X 10"* 


NGC4342 


SO 


11.4 


225 


2.2(1.3,0.8) 


X 10* 


s-23 


0.34 


0.8 


-16.4 


-20.1 


-20.7 


-20.7 


0.29 ±0.07 


1.0 ±0.3 X 10"* 


NGC0821 


E4 


24.1 


209 


3.7(2.4,0.8) 


X 10' 


s-4 


0.03 


0.6 


-20.4 


-24.4 


-24.9 


-24.8 


20±5 


6.2 ± 1.7 X 10" 


NGC2778 


E2 


22.9 


175 


1.4(0.8,0.9) 


X 10' 


s-4 


0.02 


0.4 


-18.6 


-22.0 


-22.8 


-23.0 


3.0 ±0.8 


6.5 ± 1.7 X 10'° 


NGC7052 


E4 


71.4 


266 


4.0(2.8,1.6) 


X 10* 


g-24 


0.07 


0.5 


-21.7 


-25.2 


-25.9 


-26.1 


12±3 


6.0 ± 1.6 X 10" 



Note. — (1) Galaxy Name. (2) Morphological type from RC3. (3) Galaxy Distance in Mpc. (4) Stellar velocity dispersions from T{)2 and Kormendy & Gebhardt 
2001 in units of km s"'. All cr, values have ±5% errors except for the Milky Way (±20 km s"'), Cygnus A (±90 km s"', Thornton, Stockton, & Ridgway 1999) and 
NGC5252 (±27 km s"' ). (5) BH mass in units of Mq . '(-f.-)' indicate the ± errors. (6) Method of Mbh determination (g-gas kinematics, s-stellar dynamics, m-gas 
kinematics with maser spots) and references from where Mbh was obtained (if necessary it was rescaled to the distances in column 3). (7) Black hole sphere of influence, 
Sbh = GMbh/ct^, in arcseconds. i^)N,es = IRBu/Xres where Jf^j is the spatial resolution of the observations (9) Absolute bulge B luminosity from T02 or extracted from 
the RC3 catalog (de Vaucouleurs et al. 1991). (10,11,12) Absolute J, H and K bulge magnitudes. Milky Way J and K values are taken from Dwek et al. (1995) while for 
M31 we have corrected the total galaxy magnitudes Malhotra et al. (1996) with the Simien & de Vaucouleurs (1986) disk-bulge average ratios: H values are interpolated. 
All M], Mh and Mk values have errors ±0.1 mag except for the Milky Way (AMj = 0.75, AMh = 0.75, AMk = 0.75), M31 (0.4, 0.5. OA), Centaurus A (0.2, 0.3, 0.2), 
NGC 1068 (0.6, 0.6, 0.6) and NGC 7457 (0.6, 0.4, 0.4). (13) J band effective bulge radius in kpc. For the Milky Way we have taken the estimate by T02, while for M31 
we have used that by Kormendy & Bender (1999). (14) Virial bulge mass in units of Mq . 

* Though the Milky Way represents by far the best case for a BH, it has been placed in group 2 because existing measurements of the bulge near-lR luminosity are 
uncertain and because it is beyond the scope of this paper to estimate the luminosity of the Milky Way bulge from 2MASS data. This is also the case for M3 1 . NGC 1068 
is in group 2 because the BH mass estimate is not 'secure' in the sense that the maser spots are moving sub-keplerianly (Greenhill et al. 1996) and Mbh depends on the 
adopted disk model (Lodato & Bertin 2003). Also the complex morphology did not allow to obtain an accurate estimate of the bulge luminosity. In the cases of NGC 
4459 and NGC 4596, the data do not allow a tight consti^nt on Mbh (Sarzi et al. 2001). 



References. — 1) Hermstein et al. 1999; 2) Macchetto et al. 1997; 3) Emsellem, Dejonghe, & Bacon 1999; 4) Gebhardt et al. 2003; 5) Devereux et al. 2003; 6) average 
of Bower et al. 1998; Maciejewski & Binney 2001; 7) Verolme et al. 2002; 8) Marconi et al. 2001; 9) average of Verdoes Klein et al. 2000 and Cappellari et al. 2002; 10) 
Capetti et al. 2003; 11) Sarzi et al. 2001; 12) Kormendy 1988; 13) Barth et al. 2001; 14) Tadhunter et al. 2003; 15) Ferrarese et al. 1996; 16) T02; 17) Gebhardt et al. 
2000; 18) Bower et al. 2001; 19) Ferrarese & Ford 1999; 20) Ghez et al. 2003; Schodel et al. 2002; 21) Tremaine 1995; Kormendy & Bender 1999; Bacon et al. 2001; 
22) Lodato & Bertin 2003; Greenhill et al. 1996; 23) Cretton & van den Bosch 1999; 24) van der Marel & van den Bosch 1998 



vations. Additional reasons for placing galaxies in Group 2 are 
given in Table 1 . 

3. IMAGE ANALYSIS 

We have constructed a homogeneous set of NIR images of 
the galaxies presented in Table 1 (except for the Milky Way and 



M31) by retrieving J, H, and K atlas images from the 2-Micron 
All Sky Survey (2MASS; http://www.ipac.caltech.edu/2mass). 
When a single atlas image contained only a portion of the 
galaxy, we also retrieved adjacent tiles and mosaicked the im- 
ages after subtracting the sky background and rescaling for the 
different zero points. The 2MASS images are photometrically 
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log kbui [Lq.kI log Mb„, [Mq] log R, [kpc] 

Fig. 1 . — Left (a): Mbh vs Z.K,bui for the galaxies of Group 1. The .solid line.s are obtained with the bisector linear regression algorithm by Akritas & Bershady 
(1996), while the dashed Unes are ordinary least square fits. Center (b): Mbh vs Mbui with the same notation as in the previous panel. Right (c): residuals of Mbh-c<! 
vs Re, in which we use the MBH-Ce regression by T02. 



calibrated with a typical accuracy of a few percent. More details 
can be found in Hunt & Marconi (2003; hereafter Paper II). 

We performed a 2D bulge/disk decomposition of the images 
using the program GALFIT (Peng et al. 2002) which is made 
publicly available by the authors. This code allows the fitting of 
several components with different functional shapes (e.g., gen- 
eralized exponential (Sersic) and simple exponential laws); the 
best fit parameters are determined by minimizing x^- More de- 
tails on GALFIT can be found in Peng et al. (2002). We fit 
separately the J, H and K images. Each fit was started by fit- 
ting a single Sersic component and constant background. When 
necessary (e.g., for spiral galaxies), an additional component 
(usually an exponential disk) was added. In many cases these 
initial fits left large residuals and we thus increased the number 
of components (see also Peng et al. 2002). The fits are de- 
scribed in detail in Paper II. In Table 1 we present the J, H and 
K bulge magnitudes, effective bulge radii Rg in the J band, and 
their uncertainties. The J, H and K magnitudes were corrected 
for Galactic extinction using the data by Schlegel, Finkbeiner, 
& Davis (1998). We used the J band to determine Rg because 
the images tend to be flatter, and thus the background is better 
determined. 

4. RESULTS AND DISCUSSION 

In Fig. 1 we plot, from left to right, Mbh vs /.K.buh Mbh 
vs Mbui, and the residuals of Mbh-Cc vs Rg (based on the fit 
from T02). Only Group 1 galaxies are shown. Mbui is the virial 
bulge mass given by kRg / G; if bulges behave as isothermal 
spheres, ^ = 8/3. However, comparing our virial estimates Mbui 
with those M^yn, obtained from dynamical modeling (Magor- 
rian et al. 1998; Gebhardt et al. 2003), shows that Mbui and 
Mdyn are well correlated (r = 0.88); setting k = 3 (rather than 
8/3) gives an average ratio of unity. Therefore, we have used 
k = 3 in the above formula. Considering the uncertainties of 
both mass estimates, the scatter of the ratio Mbui/Mjyn is 0.21 
dex. We fit the data with the bisector Unear regression from 
Akritas & Bershady (1996) which allows for uncertainties on 
both variables and intrinsic dispersion. The FITEXY routine 
(Press et al. 1992) used by T02 gives consistent results (see 
Fig. 1). Fit results of Mbh vs galaxy properties for Group 1 and 
the combined samples are summarized in Table 2. The intrin- 
sic dispersion of the residuals (rms) has been estimated with 



a maximum likelihood method assuming normally-distributed 
values. Inspection of Fig. 1 and Table 2 show that LK,bui and 
Mbui correlate well with the BH mass. The correlation between 
Mbh and Mbui is equivalent to that between the radius of the BH 
sphere of influence Rbh (= GMbh/ crj^ and Rg. 

4.1. Intrinsic Dispersion of the Correlations 

To compare the scatter of MBH-ibui for different wavebands, 
we have also analyzed the B-band bulge luminosities for our 
sample. The upper limit of the intrinsic dispersion of the Mbh- 
Lbui correlations goes from ~ 0.5 dex in logMBH when con- 
sidering all galaxies, to ~ 0.3 when considering only those of 
Group 1. Hence, for galaxies with reliable Mbh and Lbui, the 
scatter of MnH-^-bui correlations is ~ 0.3, independently of the 
spectral band used (B or JHK), comparable to that of Mbh-Cc- 
This scatter would be smaller if measurement errors are under- 
estimated. McLure & Dunlop (2002) and Erwin et al. (2003) 
reached a similar conclusion using R-band Lbui but on smaller 
samples. The correlation between R-band bulge light concen- 
tration and Mbh has a comparable scatter (Graham et al. 2001). 

Since MBH-ie.bui and MBH-^NiR.bui have comparable disper- 
sions, the rough bulge/disk decomposition (§1), the larger red- 
dening and stellar population effects do not apparently compro- 
mise the correlation. Most of the galaxies in the sample are 
early types, and thus may be less sensitive to the above effects. 
However, the scatter in the MBH-^bui correlations does not de- 
crease significantly when considering only elliptical galaxies. 

The correlation between Mbh and Mbui has a slightly lower 
dispersion (0.25 versus 0.3) than MBH-i'bui- If the scatter of 
Mbui-Mdyn (0.21 dex) is an indication of the additional uncer- 
tainties on our virial estimates, then the intrinsic scatter of Mbh- 
Mbui drops to ~ 0.15 dex. Judging from the present data where 
only secure Mbh and Lbui are included, ag, Mbui, Z^B.bui and 
^NiR.bui provide equally good Mbh estimates to within a factor 
of ~2. 

4.2. Correlation slopes 

All the slopes are roughly unity, but those of MBH-^bui are 
systematically larger than that of MBH-Mbui- This is expected 
if Lbui correlates with Mbh because of its dependence on Mbui 
through the stellar M/L ratio. From our Mbui-Z^bui relation, we 
find that log Tk = 0.18 logLK,bui - 2 . 1 ; the weak dependence of 
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Table 2 

Fl l' RHSULI S (logMtiH =(l + hX). 



X 




Group 1 Galaxies 




AH Galaxies 






a 


b 


rms 


a 


b 


rms 


logi-B.bur 


-10.0 


8.18 ±0.08 


1.19±0.12 


0.32 


8.07 ±0.09 


1.26±0.13 


0.48 


logij,bui- 


-10.7 


8.26 ±0.07 


1.14±0.12 


0.33 


8.10±0.10 


1.24 ±0.15 


0.53 


logi-H.bul 


-10.8 


8.19 ±0.07 


1.16±0.12 


0.33 


8.04±0.10 


1.25 ±0.15 


0.52 


logi-K.bul 


-10.9 


8.21 ±0.07 


1.13±0.12 


0.31 


8.08 ±0.10 


1.21±0.13 


0.51 


logMbul- 


10.9 


8.28 ±0.06 


0.96 ±0.07 


0.25 


8. 12 ±0.09 


1.06 ±0.10 


0.49 



T on L fully accounts for the different slopes of MeH-i'K.bui and 
M-BH-Mhu\, and the same applies to the J and H bands. 

All correlations are thus consistent with a direct proportion- 
ality between Mbh and bulge mass. This contrasts with pre- 
vious claims of a non-linearity of the MBH-Mbu\ relation (Laor 
2001) but is in agreement with McLure & Dunlop (2002). A 
partial correlation analysis of logMen (variable Xi), logCe {x2), 
and logRe (xa) shows that Mbh is separately significantly cor- 
related both with (Te and The Pearson partial correlation 
coefficients, in which the known dependence of a-g and is 
eliminated, are rn = 0.83, ris = 0.65, with a significance of 
> 99.9%. This is shown graphically in Fig. Ic, where the resid- 
uals of the T02 Msu-ag correlation are plotted against Rg ', there 
is a weak, but significant, correlation of these residuals with Rg. 
Consequently, when galaxy structural parameters are measured 
carefully from 2D image analysis, the additional, weaker, de- 
pendence of Mbh on Rg is uncovered. Thus, a combination of 
both cTg and Rg is necessary to drive the correlations between 
Mbh and other bulge properties. This fundamental plane of 



black holes will be further investigated elsewhere. 

The average logMBn/Mbui can be estimated assuming a log- 
normal distribution with normally distributed observational er- 
rors. With maximum Ukelihood we find (logMBH/Mbui) = 
-2.63 with an intrinsic dispersion of 0.27 dex (-2.79 and 0.49 
dex for all galaxies). Adopting the method of Merritt & Fer- 
rarese (2001), we find (logMBH/Mbui) = -2.81 and rms = 0.36 
(-2.86 and 0.44 for all galaxies) consistently with their result of 
-2.9 and 0.45 dex (see also McLure & Dunlop 2002). 
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